Narrow transmission bottlenecks and limited within-host viral diversity during a SARS-CoV-2 outbreak on a fishing boat (preprint)

The long-term evolution of viruses is ultimately due to viral mutants that arise within infected individuals and transmit to other individuals. Here we use deep sequencing to investigate the transmission of viral genetic variation among individuals during a SARS-CoV-2 outbreak that infected the vast majority of crew members on a fishing boat. We deep-sequenced nasal swabs to characterize the within-host viral population of infected crew members, using experimental duplicates and strict computational filters to ensure accurate variant calling. We find that within-host viral diversity is low in infected crew members. The mutations that did fix in some crew members during the outbreak are not observed at detectable frequencies in any of the sampled crew members in which they are not fixed, suggesting viral evolution involves occasional fixation of low-frequency mutations during transmission rather than persistent maintenance of within-host viral diversity. Overall, our results show that strong transmission bottlenecks dominate viral evolution even during a superspreading event with a very high attack rate.

Clinical performance characteristics of the Swift Normalase Amplicon Panel for sensitive recovery of SARS-CoV-2 genomes (preprint)

Amplicon-based sequencing methods have been central in characterizing the diversity, transmission and evolution of SARS-CoV-2, but need to be rigorously assessed for clinical utility. Here, we validated the Swift Biosciences SARS-CoV-2 Swift Normalase Amplicon Panels using remnant clinical specimens. High quality genomes meeting our established library and sequence quality criteria were recovered from positive specimens with a 95% limit of detection of 40.08 SARS-CoV-2 copies/PCR reaction. Breadth of genome recovery was evaluated across a range of Ct values (11.3 - 36.7, median 21.6). Out of 428 positive samples, 406 (94.9%) generated genomes with < 10% Ns, with a mean genome coverage of 13,545X/SD 8,382X. No genomes were recovered from PCR-negative specimens (n = 30), or from specimens positive for non-SARS-CoV-2 respiratory viruses (n = 20). Compared to whole-genome shotgun metagenomic sequencing (n = 14) or Sanger sequencing for the spike gene (n = 11), pairwise identity between consensus sequences was 100% in all cases, with highly concordant allele frequencies (R2 = 0.99) between Swift and shotgun libraries. When samples from different clades were mixed at varying ratios, expected variants were detected even in 1:99 mixtures. When deployed as a clinical test, 268 tests were performed in the first 23 weeks with a median turnaround time of 11 days, ordered primarily for outbreak investigations and infection control.

Host-pathogen dynamics in longitudinal clinical specimens from patients with COVID-19 (preprint)

Rapid dissemination of SARS-CoV-2 sequencing data to public repositories has enabled widespread study of viral genomes, but studies of longitudinal specimens from infected persons are relatively limited. Analysis of longitudinal specimens enables understanding of how host immune pressures drive viral evolution in vivo. Here we performed sequencing of 49 longitudinal SARS-CoV-2-positive samples from 20 patients in Washington State collected between March and September of 2020. Viral loads declined over time with an average increase in RT-PCR cycle threshold (Ct) of 0.87 per day. We found that there was negligible change in SARS-CoV-2 consensus sequences over time, but identified a number of nonsynonymous variants at low frequencies across the genome. We observed enrichment for a relatively small number of these variants, all of which are now seen in consensus genomes across the globe at low prevalence. In one patient, we saw rapid emergence of various low-level deletion variants at the N-terminal domain of the spike glycoprotein, some of which have previously been shown to be associated with reduced neutralization potency from sera. In a subset of samples that were sequenced using metagenomic methods, differential gene expression analysis showed a downregulation of cytoskeletal genes that was consistent with a loss of ciliated epithelium during infection and recovery. We also identified co-occurrence of bacterial species in samples from multiple hospitalized individuals. These results demonstrate that the intrahost genetic composition of SARS-CoV-2 is dynamic during the course of COVID-19, and highlight the need for continued surveillance and deep sequencing of minor variants.

Phylogenetic estimates of SARS-CoV-2 introductions into Washington State (preprint)

Background: The first confirmed case of SARS-CoV-2 in North America was identified in Washington state on January 21, 2020. We aimed to quantify the number and temporal trends of out-of-state introductions of SARS-CoV-2 into Washington. Methods: We conducted a phylogenetic analysis of 11,422 publicly available whole genome SARS-CoV-2 sequences from GISAID sampled between December 2019 and September 2020. We used maximum parsimony ancestral state reconstruction methods on time-calibrated phylogenies to enumerate introductions/exports, their likely geographic source (e.g. US, non-US, and between eastern and western Washington), and estimated date of introduction. To incorporate phylogenetic uncertainty into our estimates, we conducted 5,000 replicate analyses by generating 25 random time-stratified samples of non-Washington reference sequences, 20 random polytomy resolutions, and 10 random resolutions of the reconstructed ancestral state. Results: We estimated a minimum 287 separate introductions (median, range 244-320) into Washington and 204 exported lineages (range 188-227) of SARS-CoV-2 out of Washington. Introductions began in mid-January and peaked on March 29, 2020. Lineages with the Spike D614G variant accounted for the majority (88%) of introductions. Overall, 61% (range 55-65%) of introductions into Washington likely originated from a source elsewhere within the US, while the remaining 39% (range 35-45%) likely originated from outside of the US. Intra-state transmission accounted for 65% and 28% of introductions into eastern and western Washington, respectively. Conclusions: There is phylogenetic evidence that the SARS-CoV-2 epidemic in Washington is continually seeded by a large number of introductions, and that there was significant inter- and intra-state transmission. Due to incomplete sampling our data underestimate the true number of introductions.

Specific allelic discrimination of N501Y and other SARS-CoV-2 mutations by ddPCR detects B.1.1.7 lineage in Washington State (preprint)

Real-time epidemiological tracking of variants of interest can help limit the spread of more contagious forms of SARS-CoV-2, such as those containing the N501Y mutation. Typically, genetic sequencing is required to be able to track variants of interest in real-time. However, sequencing can take time and may not be accessible in all laboratories. Genotyping by RT-ddPCR offers an alternative to sequencing to rapidly detect variants of concern through discrimination of specific mutations such as N501Y that is associated with increased transmissibility. Here we describe the first cases of the B.1.1.7 lineage of SARS-CoV-2 detected in Washington State by using a combination of RT-PCR, RT-ddPCR, and next-generation sequencing. We screened 1,035 samples positive for SARS-CoV-2 by our CDC-based laboratory developed assay using ThermoFishers multiplex RT-PCR COVID-19 assay over four weeks from late December 2020 to early January 2021. S gene dropout candidates were subsequently assayed by RT-ddPCR to confirm four mutations within the S gene associated with the B.1.1.7 lineage: a deletion at amino acid (AA) 69-70 (ACATGT), deletion at AA 145, (TTA), N501Y mutation (TAT), and S982A mutation (GCA). All four targets were detected in two specimens, and follow-up sequencing revealed a total of 10 mutations in the S gene and phylogenetic clustering within the B.1.1.7 lineage. As variants of concern become increasingly prevalent, molecular diagnostic tools like RT-ddPCR can be utilized to quickly, accurately, and sensitively distinguish more contagious lineages of SARS-CoV-2.

Viral genomes reveal patterns of the SARS-CoV-2 outbreak in Washington State (preprint)

The rapid spread of SARS-CoV-2 has gravely impacted societies around the world. Outbreaks in different parts of the globe are shaped by repeated introductions of new lineages and subsequent local transmission of those lineages. Here, we sequenced 3940 SARS-CoV-2 viral genomes from Washington State to characterize how the spread of SARS-CoV-2 in Washington State (USA) was shaped by differences in timing of mitigation strategies across counties, as well as by repeated introductions of viral lineages into the state. Additionally, we show that the increase in frequency of a potentially more transmissible viral variant (614G) over time can potentially be explained by regional mobility differences and multiple introductions of 614G, but not the other variant (614D) into the state. At an individual level, we see evidence of higher viral loads in patients infected with the 614G variant. However, using clinical records data, we do not find any evidence that the 614G variant impacts clinical severity or patient outcomes. Overall, this suggests that at least to date, the behavior of individuals has been more important in shaping the course of the pandemic than changes in the virus.

SARS-CoV-2 ORF6 disrupts nucleocytoplasmic transport through interactions with Rae1 and Nup98 (preprint)

RNA viruses that replicate in the cytoplasm often disrupt nucleocytoplasmic transport to preferentially translate their own transcripts and prevent host antiviral responses. The Sarbecovirus accessory protein ORF6 has previously been shown to be the major inhibitor of interferon production in both SARS-CoV and SARS-CoV-2. SARS-CoV-2 ORF6 was recently shown to co-purify with the host mRNA export factors Rae1 and Nup98. Here, we demonstrate SARS-CoV-2 ORF6 strongly represses protein expression of co-transfected reporter constructs and imprisons host mRNA in the nucleus, which is associated with its ability to co-purify with Rae1 and Nup98. These protein-protein interactions map to the C-terminus of ORF6 and can be abolished by a single amino acid mutation in Met58. Overexpression of Rae1 restores reporter expression in the presence of SARS-CoV-2 ORF6. We further identify an ORF6 mutant containing a 9-amino acid deletion, ORF6 {Delta}22-30, in multiple SARS-CoV-2 clinical isolates that can still downregulate the expression of a co-transfected reporter and interact with Rae1 and Nup98. SARS-CoV ORF6 also interacts with Rae1 and Nup98. However, SARS-CoV-2 ORF6 more strongly co-purifies with Rae1 and Nup98 and results in significantly reduced expression of reporter proteins compared to SARS-CoV ORF6, a potential mechanism for the delayed symptom onset and pre-symptomatic transmission uniquely associated with the SARS-CoV-2 pandemic. ImportanceSARS-CoV-2, the causative agent of COVID-19, is an RNA virus with a large genome that encodes accessory proteins. While these accessory proteins are not required for growth in vitro, they can contribute to the pathogenicity of the virus. One of SARS-CoV-2s accessory proteins, ORF6, was recently shown to co-purify with two host proteins, Rae1 and Nup98, involved in mRNA nuclear export. We demonstrate SARS-CoV-2 ORF6 interaction with these proteins is associated with reduced expression of a reporter protein and accumulation of poly-A mRNA within the nucleus. SARS-CoV ORF6 also shows the same interactions with Rae1 and Nup98. However, SARS-CoV-2 ORF6 more strongly represses reporter expression and co-purifies with Rae1 and Nup98 compared to SARS-CoV ORF6. The ability of SARS-CoV-2 ORF6 to more strongly disrupt nucleocytoplasmic transport than SARS-CoV ORF6 may partially explain critical differences in clinical presentation between the two viruses.

In vivo antiviral host response to SARS-CoV-2 by viral load, sex, and age (preprint)

Despite limited genomic diversity, SARS-CoV-2 has shown a wide range of clinical manifestations in different patient populations. The mechanisms behind these host differences are still unclear. Here, we examined host response gene expression across infection status, viral load, age, and sex among shotgun RNA-sequencing profiles of nasopharyngeal swabs from 430 individuals with PCR-confirmed SARS-CoV-2 and 54 negative controls. SARS-CoV-2 induced a strong antiviral response with upregulation of antiviral factors such as OAS1-3 and IFIT1-3, and Th1 chemokines CXCL9/10/11, as well as a reduction in transcription of ribosomal proteins. SARS-CoV-2 culture in human airway epithelial cultures replicated the in vivo antiviral host response. Patient-matched longitudinal specimens (mean elapsed time = 6.3 days) demonstrated reduction in interferon-induced transcription, recovery of transcription of ribosomal proteins, and initiation of wound healing and humoral immune responses. Expression of interferon-responsive genes, including ACE2, increased as a function of viral load, while transcripts for B cell-specific proteins and neutrophil chemokines were elevated in patients with lower viral load. Older individuals had reduced expression of Th1 chemokines CXCL9/10/11 and their cognate receptor, CXCR3, as well as CD8A and granzyme B, suggesting deficiencies in trafficking and/or function of cytotoxic T cells and natural killer (NK) cells. Relative to females, males had reduced B and NK cell-specific transcripts and an increase in inhibitors of NF-{kappa}B signaling, possibly inappropriately throttling antiviral responses. Collectively, our data demonstrate that host responses to SARS-CoV-2 are dependent on viral load and infection time course, with observed differences due to age and sex that may contribute to disease severity.

Identification of multiple large deletions in ORF7a resulting in in-frame gene fusions in clinical SARS-CoV-2 isolates (preprint)

Peculiar among human RNA viruses, coronaviruses have large genomes containing accessory genes that are not required for replication. Numerous mutations within the SARS-CoV-2 genome have been described but few deletions in the accessory genes of SARS-CoV-2 have been reported. Here, we report two large deletions in ORF7a, both of which produce new open reading frames (ORFs) through the fusion of the N-terminus of ORF7a and a downstream ORF.

Cryptic transmission of SARS-CoV-2 in Washington State (preprint)

Following its emergence in Wuhan, China, in late November or early December 2019, the SARS-CoV-2 virus has rapidly spread throughout the world. On March 11, 2020, the World Health Organization declared Coronavirus Disease 2019 (COVID-19) a pandemic. Genome sequencing of SARS-CoV-2 strains allows for the reconstruction of transmission history connecting these infections. Here, we analyze 346 SARS-CoV-2 genomes from samples collected between 20 February and 15 March 2020 from infected patients in Washington State, USA. We found that the large majority of SARS-CoV-2 infections sampled during this time frame appeared to have derived from a single introduction event into the state in late January or early February 2020 and subsequent local spread, strongly suggesting cryptic spread of COVID-19 during the months of January and February 2020, before active community surveillance was implemented. We estimate a common ancestor of this outbreak clade as occurring between 18 January and 9 February 2020. From genomic data, we estimate an exponential doubling between 2.4 and 5.1 days. These results highlight the need for large-scale community surveillance for SARS-CoV-2 introductions and spread and the power of pathogen genomics to inform epidemiological understanding.

Coast-to-coast spread of SARS-CoV-2 in the United States revealed by genomic epidemiology (preprint)

Since its emergence and detection in Wuhan, China in late 2019, the novel coronavirus SARS-CoV-2 has spread to nearly every country around the world, resulting in hundreds of thousands of infections to date. The virus was first detected in the Pacific Northwest region of the United States in January, 2020, with subsequent COVID-19 outbreaks detected in all 50 states by early March. To uncover the sources of SARS-CoV-2 introductions and patterns of spread within the U.S., we sequenced nine viral genomes from early reported COVID-19 patients in Connecticut. Our phylogenetic analysis places the majority of these genomes with viruses sequenced from Washington state. By coupling our genomic data with domestic and international travel patterns, we show that early SARS-CoV-2 transmission in Connecticut was likely driven by domestic introductions. Moreover, the risk of domestic importation to Connecticut exceeded that of international importation by mid-March regardless of our estimated impacts of federal travel restrictions. This study provides evidence for widespread, sustained transmission of SARS-CoV-2 within the U.S. and highlights the critical need for local surveillance.

Comparative Performance of SARS-CoV-2 Detection Assays using Seven Different Primer/Probe Sets and One Assay Kit (preprint)

More than 100,000 people worldwide are known to have been infected with SARS-CoV-2 beginning in December 2019. The virus has now spread to over 93 countries including the United States, with the largest cluster of US cases to date in the Seattle metropolitan area in Washington. Given the rapid increase in the number of local cases, the availability of accurate, high-throughput SARS-CoV-2 testing is vital to efforts to manage the current public health crisis. In the course of optimizing SARS-CoV-2 testing performed by the University of Washington Clinical Virology Lab (UW Virology Lab), we tested assays using seven different primer/probe sets and one assay kit. We found that the most sensitive assays were those the used the E-gene primer/probe set described by Corman et al. (Eurosurveillance 25(3), 2020, https://doi.org/10.2807/1560-7917.ES.2020.25.3.2000045) and the N2 set described by the CDC (Division of Viral Diseases, Centers for Disease Control and Prevention, 2020, https://www.cdc.gov/coronavirus/2019-ncov/downloads/rt-pcr-panel-primer-probes.pdf). All assays tested were found to be highly specific for SARS-CoV-2, with no cross-reactivity with other respiratory viruses observed in our analyses regardless of the primer/probe set or kit used. These results will provide invaluable information to other clinical laboratories who are actively developing SARS-CoV-2 testing protocols at a time when increased testing capacity is urgently needed worldwide.